Making vortices in dipolar spinor condensates via rapid adiabatic passage 
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We propose to the create vortices in spin-1 condensates via magnetic dipole-dipole interaction. 
Starting with a polarized condensate prepared under large axial magnetic held, we show that by 
gradually inverting the held, population transfer among different spin states can be realized in a 
controlled manner. Under optimal condition, we generate a doubly quantized vortex state containing 
nearly all atoms in the condensate. The resulting vortex state is a direct manifestation of the dipole- 
dipole interaction and spin textures in spinor condensates. We also point out that the whole process 
can be qualitatively described by a simple rapid adiabatic passage model. 
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Ever since the first realization of vortices in an atomic 
condensate by dynamical phase imprinting [l[ , quantized 
vortices in quantum gases have attracted great attention. 
A variety of techniques have been used to generate vor- 
tices, including mechanically stirring the atomic cloud 
with laser beams Q, rotating asymmetric traps slic- 
ing through the condensate with a fast-moving perturba- 
tion 0, topological phase imprinting Q, decay of soli- 
tons [5j , coherently transferring orbital angular momen- 
tum from photons to atoms [7|, and merging multiply 
trapped condensates Q. None of these schemes, however, 
relies on the specific forms of atom-atom interactions. 

In the present work, we propose a novel scheme to 
create vortices in a spin-1 condensate by utilizing the 
magnetic dipole-dipole interaction such that the result- 
ing vortices become a direct manifestation of the under- 
lying dipolar interaction. For convenience, the three spin 
components of a spin-1 atom are labeled as a — and 
±1. Starting with a pure condensate of a = — 1 atoms 
prepared under an axial magnetic field along negative z- 
axis, we show that a doubly quantized vortex state in 
a = 1 component can be created by gradually invert- 
ing the magnetic field from negative to positive. With 
a careful control of the sweeping rate of the magnetic 
field and other parameters, the efficiency of atom trans- 
fer from a = —1 to 1 component can approach unity. 
Quite remarkably, as we will show, this dynamical evo- 
lution can be understood as a rapid adiabatic passage 
process described by a simple Landau-Zener tunneling 
model. Furthermore, due to their very rich physical prop- 
erties d H3, HI E3, EI H d, dipolar spinor conden- 
sates have become one of the focuses in the study of quan- 
tum gases. In the experiment performed at Stamper- 
Kurn's group [l^ . evidences supporting dipolar effects 
show up in the form of intriguing spin textures whose 
detection, however, requires sophisticated imaging tech- 
niques. While in our scheme, the vortex state is a much 
more robust signal and can be readily detected via a 
straightforward density measurement. 

We consider a condensate of N spin-1 87 Rb atoms. 
Within the mean- field treatment, the dynamical behavior 



of the condensate wave functions ip a (r) is described by 
(here and henceforth, summation over repeated indices 
is assumed) 
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[T + U + c p(r)]Va + gAisBcff (r) ■ F a/S V/3, (1) 



where T = — fi, 2 V 2 /(2m) with m being the mass of the 
atom, U(r) = |mw§(i 2 + y 2 + X 2 z 2 ) is the trapping po- 
tential with A being the trap aspect ratio, p(r) = ip^a is 
the total density of the condensate, and cq = 47rft 2 (ao + 
202) /(3m) characterizes the spin independent collisional 
interaction with a/ (/ = 0,2) being the s-wave scatter- 
ing length in the combined symmetric channel of total 
spin /. For rubidium atoms, we have ao = 101. 8as and 
02 = 100. Aolb with as being the Bohr radius. Further- 
more, g(— —1/2) is the Lande g-factor, \ib is the Bohr 
magneton, and F is the angular momentum operator. 
The effective field includes the external magnetic field 
B = B(t)z, and the mean fields originating from the 
spin-exchange and dipole-dipole interactions 
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where c 2 = 47rft 2 (a 2 - a )/(3m), S(r) = ip%F a pi/)p, e = 
(r — r')/|r — r'|, and the strength of dipolar interaction 
is characterized by Cd = ^o^g 2 / (4tt) with being the 
vacuum magnetic permeability. 

For the numerical results presented in this work, the 
transverse trapping frequency is taken to be luq = 
(27r)100Hz. We shall focus on three different trap ge- 
ometries corresponding to prolate, spherical, and oblate 
trapping potentials with, respectively, A = 0.25, 1, and 
6. Unless otherwise stated, the total number of atoms is 
chosen to be N = 2 x 10 6 . 

Let us first recall the ground state structure in mag- 
netic field. For convenience, we assume that the magnetic 
field is along the negative z-axis, i.e., B < 0. In an oblate 
trap, as shown in Ref. [T2|, the ground state wave func- 
tion takes the form ip a ( r ) = y Pa( r ) exp[i(w Q (p + <p a )], 
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FIG. 1: Schematic plot of the magnetic field dependence of 
the ground state structure in an oblate trap. 



where the densities p a = IV'al 2 are axially symmetric, 
w a are the winding numbers, tp is the azimuthal angle, 
and <p a are phase angles satisfying 2ifo — (pi — tp-i = 0. 
There exists a critical magnetic field strength B* such 
that (wi,w ,w-i) = (-1,0,1) and (-2,-1,0) for > 
B > —B* and B < —B*, respectively. Across the criti- 
cal field — B* , the total energy is continuous, indicating a 
second-order phase transition. Obviously, for sufficiently 
strong magnetic field, the system will be polarized. We 
can define the saturation field strength B sat > such 
that when B < —B sai , over 99.9% of the population will 
be in a = — 1 component. The saturation field is an in- 
creasing function of N and A. For parameters adopted 
in this paper, 73 sat is about several tens of micro-Gauss. 
Finally, incorporating the results corresponding to B > 
case gives us the complete picture of ground state phases 
in an oblate trap, schematically plotted in Fig. [T] 

In a spherical trap, the axial symmetry of the den- 
sities p a is lost under weak magnetic field. However, 
if one increases the magnetic field strength, p a recovers 
axial symmetry, especially when the system becomes po- 
larized. The saturation field in a spherical trap is only 
about several micro-Gauss, much lower than that in an 
oblate trap. Finally, only polarized phase exists in a pro- 
late trap. 

To study the dynamic properties, we numerically 
evolve Eq. ([T]) with the initial wave functions being the 
ground state under field Bq < — _B sat . The magnetic field 
is assumed to vary linearly as 
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FIG. 2: (Color online) Typical dynamic behaviors in an oblate 
trap (A = 6). (a) The time dependence of reduced atom 
number m for various control parameters (N, Bo,vb)- (b) 
The integrated densities p\(x,y) and phases (insets) of wave 
function tpi(x, y, 0) for (N,Bo) = (2 x 10 6 ,— 40) with various 
Vb's. (c) Same as (b) except for N — 10 7 , Bo = — 100/itG, 
and vb = 0.5. 



B(t) =B + v B t, 

where vb > (in units of /iG -ujq) is the sweeping rate of 
the magnetic field. We shall explore how to control the 
dynamic behaviors of rubidium condensates by tuning 
parameters A, N, Bq, and Vb- 

In Fig. (a) , we plot the time dependence of the re- 
duced atom number n a = J drp a (r) for a = 1 spin 
component in an oblate trap. Since n\ is essentially zero 
at the beginning, its final value can be regarded as a mea- 
sure of the atom transfer efficiency. n\ remains negligible 
until B is close to zero. It then grows with oscillations 
and eventually reaches some steady-state value. Given 
N and Bq, the asymptotic transfer efficiency increases as 
one lowers the sweeping rate, and it can be as high as 
80% for TV = 2 x 10 6 and v B = 0.15. However the onset 
of population transfer occurs earlier for larger vb- As 
we shall show below, the relation between n\ and vb can 
be understood using a simple Landau-Zener tunneling 
(LZT) model of a three-level system. 



Figure [5] (b) shows the time dependence of integrated 
density 



Pa(x,y) 




\\ i see that the density depletion appears at the center of 
the condensate ipx . Further examination of the phases of 
the wave function confirms that it is a doubly quantized 
vortex which is unstable against splitting into two singly 
quantized vortices [Tj}. Indeed one can see that, as the 
system continues to evolve, the dynamical instability of 
the wave function ipi sets in and the doubly quantized 
vortex breaks into two vortices. The time that this break 
happens is sensitive to the magnetic field sweeping rate: 
It happens at a later time for a lower sweeping rate (i.e., 
smaller vb)- 

In Fig. [2] (c), we show the evolution of the inte- 
grated density p\ for an increased total number of atoms 
TV = 10 7 . In this case, the atom transfer efficiency be- 
comes nearly unity [see Fig. [5] (a)] even for a rather large 
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FIG. 3: (Color online) Time dependence of relative popula- 
tion m (a) and integrated density pi (b) in various trapping 
potentials for Bo = — 60^iG and Vb = 0.5. Insets in (b) show 
the phases of wave function ipi(x, y, 0). 

sweeping rate vb = 0.5 and a stronger initial magnetic 
field B = — 100 ^G. The doubly quantized vortex state 
survives for a much longer time compared with the pre- 
vious case with smaller N but the same sweeping rate. 

The creation of the doubly quantized vortex can be 
most easily understood by noting that the dipolar inter- 
action conserves the total angular momentum J = F + L. 
When an axial magnetic field is present, J z is still con- 
served. For the initial condensate under study, we have 
mj = rap = — 1. If we assume that all atoms are trans- 
ferred to a = 1 state at the end of the process, the spin 
angular momentum then becomes wif = 1. As a result, 
wave function if)\ must carry an orbital angular momen- 
tum mj, = —2, representing a doubly quantized vortex, 
in order to conserve J z . 

Next, we turn to study the effects of trap geometry on 
the population transfer and the wave functions. In Fig. [31 
we present the simulation results for all three trap geome- 
tries with Bq = — 60/iG and vb — 0.5. Clearly, for given 
Bq and vb, n\ increases with trap aspect ratio A. From 
the plots of the density profile in Fig. [3] (b), it can be 
seen that the doubly quantized vortex in a prolate trap 
is most stable against splitting. In such prolate trap, the 
population transfer efficiency is low. The vortex core is 
filled by a large number of spin a = — 1 atoms, which pro- 
vides an effective pinning effect that helps to stablize the 
vortex state [l8l |. However, in both prolate and spherical 



traps, the vortices suffer a different kind of instability: 
The vortex line tends to be distorted and cannot remain 
straight along the axial direction. A consequence of such 
instability is to reduce the contrast of the vortex core 
in the integrated density profile p(x,y). From these re- 
sults, we conclude that an oblate trap provides the best 
candidate to study this phenomenon. 

To gain more insights into the dynamical properties of 
the system, let us consider the spin dynamics of a single 
spin-1 particle under external magnetic field. For the 
simplest case, it reduces to the model Hamiltonian 

#lzt = b(t)F z + b x F x , (3) 

where, without loss of generality, the transverse field b x 
is assumed to be a constant and along the x-axis. Here 
all quantities are taken to be dimensionless. The energy 
spectrum of Hlzt is schematically plotted in Fig. Q] (a) 
as a function of the axial magnetic field strength b. The 
presence of b x changes the three-level crossing at b = 
into anti-crossings. The minimum gap between two 
eigenvalues of -Hlzt is A = \/2b x . Equation ([SJ describes 
essentially the LZT of three levels . Figure H (b) 
illustrates an example of the dynamics according to the 
LZT model. Under a large initial axial field bo < 0, the 
atom is prepared in a = — 1 state. As we sweep the axial 
field linearly, bit) = bo+Vbt, over the anti-crossings, a = 
state is first populated, followed by a = 1 state. Namely, 
the wave function bifurcates into the states forming the 
anti-crossing and thus becomes a coherent superposition 
of different spin states. This coherent superposition is the 
cause for the population oscillations which are suppressed 
at large t limit. The asymptotical transition probability 
to a = 1 state is [l9[ 

Pl = (l-e-^^) 2 , (4) 

an increasing function of b x while a decreasing one of 
Vb- Moreover, the population oscillations disappear com- 
pletely for sufficiently small sweeping rate Vb as the adia- 
batic limit is reached, realizing a rapid adiabatic passage 
process which represents an important method of trans- 
ferring population from one quantum state to another. 

In Fig. [4] (c), we present the B dependence of the re- 
duced atom number n a in a dipolar spinor condensate. 
Here the transverse field, B^, is the transverse compo- 
nent of the effective field. Consequently, it has non-trivial 
spatial and temporal dependences. Despite the complex- 
ity of the condensate system, we still see a remarkable 
qualitative agreement between the full model and the 
LZT model: the latter captures all the essential features 
of the condensate population dynamics. For the parame- 
ters used in our full numerical simulations, we always see 
oscillations in the population dynamics no matter how 
slowly the magnetic field is swept, meaning we cannot 
reach the complete adiabatic limit in our full simulation. 
This is probably due to the accumulated numerical noise. 
In practice, the lifetime of the condensate (typically on 
the order of tens of seconds) may also set a constraint 
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FIG. 4: (Color online) (a) Eigenenergies of J/lzt as a function 
of axial field b. (b) Population dynamics of the LZT model 
with b x = 1, bo = —10, and Vb = 1.5. (c) Reduced atom 
numbers n a and total energy per atom £ versus magnetic field 
in an oblate trap (A = 6) . Other parameters are Bq — —40 jiG 
and vb = 0.25. 

on whether adiabatic limit can be reached. Note that 
in our simulation, we always limit the total time to be 
about 1 second. In Fig. Q] (c), we also plot the total en- 
ergy per atom £ as a function of B. Again we see that 
this curve agrees qualitatively with the lowest adiabatic 
energy level in the LZT model as shown in Fig. 2] (a). 

The transition probability Eq. ^ can also be used 
to interpret the trap geometry dependence of popula- 
tion transfer efficiency ri\. Even through \ci\ is about 
ten times larger than Cd for rubidium atom, we find in 



our simulation that the contribution to H^s [see Eq. ([2}] 
is mainly provided by the dipole-dipole interaction, and 
is enhanced by the oblate geometry. Furthermore, the 
dipolar interaction can also be enhanced by increasing 
the total number of atoms N, which explains the atom 
number dependence of the transfer efficiency n±. 

We remark that in previous studies on the Einstcin-dc 
Haas effect of dipolar spinor condensates, the magnetic 
field is usually inverted suddenly [IH, [3, [l5| . While in 
this work, the inverting process of the magnetic field is 
controlled deliberately, such that we may take advantage 
of the LZT to realize much higher population transfer 
efficiency and better control of the system. 

In conclusion, we have studied the dynamics of a ru- 
bidium spinor condensate under a time-dependent axial 
field. We show that by sweeping the magnetic field and 
inverting its direction, a doubly quantized vortex con- 
taining nearly all atoms in the condensates is created. 
Despite of the complicated nature of the full system, the 
population dynamics discussed here can be understood 
as a rapid adiabatic passage process described by a sim- 
ple three-level Landau-Zener tunneling model. Since our 
scheme critically depends on the properties of the dipo- 
lar interaction, not only does it provide a new method 
for generating vortices in atomic condensates, it can also 
be used as a simple and robust mechanism for the exper- 
imental demonstration of dipolar effects in spinor con- 
densates. We hope our work will stimulate experimental 
efforts along this line. Our future works will also include 
the study of the controlled spin dynamics of the spin-3 
chromium condensate. 
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